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Abstract 
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scaling of the renormalized variables we obtain, as fixed points of the 
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interest of a renormalization group approach to the generalized central 
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1 Introduction 



The classical Central Limit Theorem states that the properly rescaled, cen- 
tered sum of n independent and identically distributed random variables with 
finite mean and variance converges to the Gaussian distribution when n — > oo. 
This universality explains the ubiquity of the Gaussian distribution in the 
applications of probability theory to many branches of science. In Physics, 
for example, the Central Limit Theorem is behind such fundamental results 
like the Maxwell-Boltzmann distribution and the microscopic interpretation 
of diffusion by Einstein [T] . 

Diffusion equations are obtained as hydrodynamic limits of continuous time 
random walks [2] where the probability distribution of the waiting-time is 
Markovian with mean waiting-time r, and the probability distribution of the 
step-size has finite variance, a. Therefore, the microscopic transport mecha- 
nism has finite characteristic length and time scales and the diffusion coeffi- 
cient is proportional to <t 2 /t. However, in the last years and in connection 
with the study of complex systems, it has been discovered that many pro- 
cesses in Physics, Biology, Economy and the Social Sciences exhibit scale-free 
transport [HI HI IS]) i- e - transport in which characteristic spatial and/or tem- 
poral scales are lacking. The so-called anomalous diffusion of these systems is 
understood in terms of non-Gaussian statistics of the underlying microscopic 
processes and modeled by means of fractional differential equations [6j [7] . The 
Generalized Central Limit Theorem [3 [9] gives all the possible limits of sums of 
(properly rescaled) independent and identically distributed random variables, 
without the hypothesis of finite variance. Precisely the limit distributions with 
infinite variance, usually called Levy distributions, are the interesting ones for 
understanding scale-free phenomena. 

On the other hand, the Renormalization Group is the main theoretical tool 
to investigate the universality that appears in different branches of Mathe- 
matics and Physics. Essentially, the Renormalization Group explains how a 
system changes when the scale of observation is modified. A Renormalization 
Group transformation usually consists in averaging certain degrees of freedom 
in a way that the original system is mapped to another with fewer degrees of 
freedom and different coupling constants. Such transformation defines a flow in 
the space of theories, with the fixed points and their linear stability properties 
giving much information about the large scale behavior of the system. 

There is an intimate relationship between the Central Limit Theorem and 
the Renormalization Group. In particular, a proof of the former can be given 
from the perspective of the latter. This is helpful in order to understand in a 
different way the mechanism of convergence of the sums of distributions, which 
are viewed as iterations of certain Renormalization Group transformations. 
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This reformulation was rigorously done by G. Jona-Lasinio [TO] for the classical 
Central Limit Theorem, see also \12\ [TO] . 

As far as we know, the aforementioned appearance of Levy statistics in a 
number of transport processes is not yet well-understood from first principles. 
Let us focus for a moment in the observed anomalous transport in certain 
regimes of turbulent fusion plasmas [14\ [TO! [TO] . The conundrum can be in- 
formally posed as follows: the fundamental equations are non-linear partial 
differential equations (fluid momentum balance equation coupled to Maxwell 
equations, for example), but particle transport seems to be suitably modeled 
in terms of linear fractional differential equations. Let us simply mention 
that whereas the symbol in Fourier space of an ordinary derivative operator 
is {—ik) n with n 6 Z, the symbol of a fractional derivative operator is {—ik) a 
with a E R. A deep understanding of this change in the analyticity properties 
of the involved operators (or equivalently, a deep understanding of the emer- 
gence of Levy statistics) is still lacking. From the point of view of Physics it is 
natural to explore the application of Renormalization Group ideas. Although 
a satisfactory and complete answer is probably far ahead, in the present paper 
we try to take a first step by showing how the Generalized Central Limit Theo- 
rem and therefore Levy distributions show up from Renormalization Group 
arguments. We put special emphasis on the study of the flow in the space 
of probability distributions. We think that this might open a way to make 
further progress in the comprehension of the problems stated above. 

We find it valuable to further motivate a Renormalization Group approach 
to the Generalized Central Limit theorem by working out rather briefly a 
simple but interesting physical model in which such a perspective emerges. 
Consider a discrete-time random walk equation 



where n(x,t) is the density of walkers and p is an arbitrary symmetric proba- 
bility density function (p.d.f. in the following), so that p(u)du is the proba- 
bility of a step taking value in [it, it + dit]. Jumps take place at discrete-time 
intervals of length r, i.e. t = nr, n £ N. The so-called fluid or hydrodynamic 
limit of this equation is obtained by considering that the previous process oc- 
curs at a microscopic scale well-separated from the macroscopic observational 
space-time scales. There are at least two ways of implementing this idea. 

The first way (see |17j for a clear and detailed treatment in more general 
random walk models) consists in relating the microscopic and macroscopic 
characteristic time and length by means of a scaling factor A that eventually 
goes to infinity. As long as the dimensions are properly chosen, (11. ip yields a 
(fractional) differential equation when A — > oo. Namely, define 




(1.1) 
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n(x, t) :=n(A A x, At), 

the density in macroscopic units, where A is the scaling factor that converts 
the microscopic scale into the macroscopic one. The dimension A has to be 
adjusted to obtain a non trivial fluid limit when A — > oo. We can write now the 
microscopic random walk equation (jl.ip in terms of the macroscopic density 

n(x,t + r/A) = A A y p(A A (x - y))n(y, t)dy. (1.2) 

Or considering its Fourier transform 

n(k,t + r/A) = / 5(A~ A k)n(k,t). (1.3) 

If we assume 

p( k ) = i-c\k\ a +o(\kn, (i.4) 

for some a 6 (0,2], c > 0, and take A = 1/a, we obtain a non-trivial limit 
when A — > oo: 

a t n(k,t) = --|k| Q n(k,t), (1.5) 

r 

that represents a diffusion equation with fractional derivatives for the density 
function n(x, t). In particular, Einstein's derivation of the diffusion equation 
from a random walk model corresponds to a = 2. 

Note that in the previous limit most of the details of the microscopic 
process have disappeared and only the lowest order in the expansion of p — 1 
around k = is relevant for determining the form of the fractional diffusion 
equation. This is a manifestation of universality, typical in phenomena that 
emerge in the macroscopic domain. Observe also that the description obtained 
at the macroscopic level, a fractional differential equation, is of different nature 
to that at the microscopic level, a discrete time random walk equation. 

The second way to tackle the problem is inspired by the Renormalization 
Group ideas. We will see that in this second approach the same description 
of the process, i.e. a discrete time random walk, is kept in the microscopic 
and macroscopic regimes. This procedure is specially useful in situations like 
quantum field theory where, due to the appearance of infinities, a cut-off is 
usually required to define the theory. To understand how this works we start 
by iterating (jl.ip to get 

n(x, t + 2r) = J p* p(x- y)n(y, t)dy, (1.6) 

where p * p(x) = J p(x — y)p{y)dy stands for the convolution of p with itself. 
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In order to compare this equation with the original one we perform a 
rescaling in space and time and introduce the new density 




n(ax', 2t') 



a > 0. 



In terms of this density Eq. (|1.6p reads 



n 




(1.7) 



that looks exactly like (jl.ip except that the p.d.f. has changed to 



T a p(x) = ap* p(ax). 



We call T a the renormalization group transformation. 

If we examine carefully what (jl.7p means, we observe that although it 
looks like a process in discrete time r with p.d.f. T a p, it actually represents 
the original process with p.d.f. p, discrete time r/2 and rescaled space variable. 
The fluid or diffusive limit is obtained when we apply this procedure infinitely 
many times, so that we can write 



Therefore, in order to have a well defined process in the fluid limit we should 
adjust a so that the previous limit exists, p is necessarily a fixed point of the 
renormalization group transformation. Universality in this approach shows up 
because the fluid limit is the same for all initial p.d.f. 's p that belong to the 
domain of attraction of the same fixed point under the action of T a . All these 
considerations motivate the interest of the study of the renormalization group 
transformation T a , its fixed points and their stability. 

The rest of the paper is organized as follows. In Section 2 we introduce 
more precisely the renormalization group transformation and study the be- 
havior of the moments. In Section 3 we establish the appropriate differential 
setup for investigating topological properties of the transformation. Finally, 
Section 4 is devoted to the study of the fixed points of the transformation, 
their stability, and their domain of attraction. 

2 Renormalization group transformation 

Let p be a p.d.f. in R, i.e. 




(1.8) 



where 



p = Urn IT p. 

m— ¥oo 
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(i) p{x) > 0, Vx G R and 



we define the following Renormalization Group transformation acting on it 



where a is a real number different from zero. Note that for normalized p, T a p 
is also normalized. 

The meaning of the transformation is clear; T a p is the probability density 
of the random variable £' = (£i+£2)/o where ^1,^2 are independent identically 
distributed random variables with density p. The two essential ingredients of 
the renormalization group: new variables that represent an average of the old 
ones and the scaling of the former, are present in this example in the simplest 
way. 

To understand how the different values of a lead to different properties 
of the transformation it is instructive to consider the situation in which all 
moments 



are finite (in the subsequent sections we will work out the general case) and 
study their behavior under the transformation. One immediately obtains the 
following expression for the moments of the transformed probability distribu- 
tion: 



which implies that the transformation maps densities with zero mean into 
themselves. In this case the relation between the variances is 



Now, we can distinguish three different regimes after m successive appli- 
cations of the transformation T a according to the values of a, when m — > 00: 

- For |a| > \/2, 





(2.1) 



i=o v ' 

If we take in particular n = 1 we get: 

2 

(x)T a p = -<s)p, 





m— ¥00 




G 



which implies that for distributions with finite variance, Tj^p should 
approach the Dirac delta function in the large m limit. In the next 
section we shall make precise the topological space in which T a acts and 
we shall see in which sense the above limit holds. 



For \a\ < y/2, assuming that (x 2 ) p / 0, we have 

lim (x 2 )T m o = oo. 



Finally, when \a\ = y/2, the value of (x 2 ) does not change under the 
application of the renormalization group transformation. In this case we 
have a fixed point with finite non-zero moments. Denoting by p the 
fixed point, T a p a = p , we will have that 



i=l v y 

which is solved by 

That solution coincides, of course, with the expression for the moments 
of the Gaussian distribution. 



3 Banach space structure 

The natural framework to ask topological and differential questions about 
T a is that of Banach spaces. In our case we shall consider test functions 
in Coo(M): the Banach space of continuous functions with vanishing limit 
at oo, endowed with the supremum norm. The probability densities are the 
positive distributions of unit norm in its topological dual C 00 (]R) / . The space 
Coo(M.y consists of the finite, complex Radon measures in M; it contains, for 
instance, distributions supported on discrete sets (Dirac delta functions) as 
well as densities in L 1 (M). 

We shall denote the action of an element of p G C <x (M.y on / G C 
pf, by 

/•oo 

Pf ■= / p(x)f(x)dx. 



It is a classical result, see [20] for instance, that the convolution p\ * p2 
of two finite Radon measures is again a finite Radon measure and actually 
C 00 (1R) / is a commutative Banach algebra with the convolution. In particular, 
the convolution is associative and satisfies 

|| Pi * P2 ||<|| Pi || || P2 ||, 
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where the equality is attained if both distributions are positive. 

To define the renormalization group transformation in C 00 (IR) / we intro- 



duce the dilation n r : Coo (R) — > Coo (R) , 

(vr r /)(x) := f(rx). 

Now we can obtain the transformed distribution under the renormalization 
group by means of the formula 

T a p = (p* p)vr 1/a (3.1) 

that correctly defines a linear functional T a p : Coo(R) — > R. 

The main differential properties of the renormalization group transforma- 
tion are collected in the following proposition. 

Proposition 1: 

The transformation T a : Coo(K)' — > C 00 (R) / in \3. 1\) is continuous and 
differentiable with bounded continuous differential. 

Proof: To show the continuity of T a it is enough to realize that %\i a is an 
isometry. Then, we get 

II T aPl - T aP2 || =|| Oi + p 2 ) * (pi - P2)n l/a || 

<|| Pi + P2 || || Pi - pi ||, (3.2) 

from which the continuity of T a follows. On the other hand, the differential of 
T a is given by 

(DT a ) p ( = 2{p * C)7r 1/a , (3.3) 
which is a bounded linear map. In fact one has 

|| (DT a ) P C || =2 || p*C||<2 || p || || C|| • 

Therefore, 

II {DT a ) p \\< 2 || p ||, 

and the equality is attained if p is positive. As for the continuity of the 
differential simply observe that 

|| (DT a ) pi - (DT a ) P2 || = || (DT a ) pi - p2 \\< 2 || pi - p 2 \\, 

and the proof is concluded. □ 
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4 Fixed points 



Much information about the renormalization group transformation is obtained 
by studying its fixed points and its linearization around them. In order to work 
out the fixed points it is convenient to express the transformation in terms of 
the characteristic functions 

fTp(k) = J T a p(x)e lkx dx = p(k/a) 2 . 

Then, for the fixed point T a p Q = p , one has 

p a (k/af = pM- (4.1) 

Note that this equation is solved by functions p whose logarithm is a symme- 
tric, homogeneous function of degree a = log 2/ log \ a\. To find more general 
solutions we distinguish two cases according to the sign of a. 

- If a is positive, we have solutions for (|4,ip of the form 

p (k) = S a , A (k) := exp(-A\k\ a 9(k) -A\k\ a 9(-k)), 

where a is given by the relation \a\ = 2 1//q , A is the complex conjugate 
of A and 9 is the Heaviside step function. 

The fact that S a ,A ls the characteristic function of a probability density 
(i.e. a positive functional) imposes additional restrictions on A and a. 
In particular it is well-known (see Refs. [El [9] for instance) that if A 7^ 0, 
then we must have 

0<a<2 and A = \A\e i,f with \tp\ < |(1 - \a - 1|). (4.2) 

The solution with A = defines a fixed point of the renormalization 
group, corresponding to the Dirac delta function, that exists for any 
value of a. 

These solutions of (|4.ip are the so-called strictly stable laws, denoted 
here by 5 QjJ 4, that where introduced by P. Levy [8], see also [9~1 fl8| [19] . 

- For negative a, the positive and negative values of k are related through 
(|4.ip . As a result, the fixed points are the ones described above for \a\ 
but restricted to real values of A. 

After discussing the fixed points of our renormalization group transforma- 
tion we analyze their stability and the domain of attraction (around the fixed 
point). The standard way to approach this problem is to study the spectrum 



and eigenvectors of the differential of the transformation at the fixed points. 
The results can be summarized in the following 

Proposition 2: 

Let (DT a ) p be the operator defined in \3. 3\) . po one of the fixed points of 
T a discussed above, and a the complex spectrum of {DT a ) po . 

(i) If p = S a< A with Re^4 > 0, then a = {A, s. t. |A| < 2} and all its 
values belong to the pointwise spectrum except for A = 0. 

(ii) If p = S a) A with Re A = 0, then a = {A, s. t. |A| = 2} with its only 
eigenvalue at A = 2. 



Proof: In order to determine the pointwise spectrum one has to solve the 
eigenvalue equation 

(DT a ) Po ( = \(. (4.3) 

We will focus first on (i), corresponding to the non-trivial fixed point po = 
S a ,A with Rej4 > 0. The equation for the eigenvalues f|4.3|) can be more easily 
handled by writing it in terms of the characteristic functions, £. Namely, 

2S aA (k/a)C(k/a) = \((k). 

Using the ansatz £(/c) = S a) A(k)rj(k) and taking into account the properties 
of the fixed point, we arrive at the following equation 

2r}(k/a) = Xr}(k), (4.4) 

which can be solved with homogeneous functions of the appropriate degree. 
Namely, for positive a = 2 1 /" we have solutions of the form 

ffc(fc) = (B + 9(k) + B_d(-k))\k\ s (4.5) 

1 s 

with eigenvalue X s = 2 «. 

In order to ensure that £ = S a ^A * V £ Coo(M)' we must have Re s > 0. We 
shall now show that it is also sufficient. The function 

U X ) = J_ r ^ k)e -A\k\"e(k)-A\k\«e(-k) e -ik X dfc> 

where Re^4 > and Res > 0, is clearly continous and even smooth, so we 
only need to show that it decays sufficiently fast for large Consider the 
integration of the positive momenta 
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F? f°° 

{ f x ) = f± / ^e-^-^dfe. 
27r Jo 

Assuming that x > and Re A > we can perform a slight clockwise 
rotation of angle 7 of the integration line in the complex plane. It does not 
affect the result due to Cauchy's theorem. If we also change to a new variable 
z = e iry kx the integral reads 

R 0— 17(3+1) poo 
Ux) = B+e ^ / z s e -A(z/a>re-^-ize-^ dz 

and for small 7 > we can safely take the limit of large x inside the integral 
to get 

C+ (x) = I®, with lim f{x) = B + ^±Ve-^y\ 

We can deal in an analogous way with the integration on the negative semiaxis 
and also with the case of negative x. In all cases we can show that provided 
Re s > the function decays fast enough so that £ = S a< A * V e CooO^)' 

In this way we prove that all values of A s.t. < |A| < 2 are in the 
pointwise spectrum. As the spectrum must be a closed set, is also in a. To 
exclude other points, i.e |A| > 2 we can use the well-known argument that if 
|A| >|| (DT a ) po ||= 2 then it belongs to the resolvent. This concludes the proof 
of part (i) of the proposition. 

We proceed to prove part (ii). There are only two possibilities of having 
a fixed point S at A with Re A = 0: either a = 1 and (p = ±ir/2 or A = 0. In 
all other instances Re ^4 > 0. We shall discuss in detail the fixed point with 
A = 0, i.e. the Dirac delta function 5, which exists for any value of a. To 
study the spectrum of the differential at this point we introduce the operator 

UC :=\{DT a ) s C = ^ 1/a . 

U is an isometry and, therefore, its spectrum lies entirely in the unit circle. 
On the other hand the equation for the eigenvalues 

= K» (4-6) 

is solved with \i = 1 and £i = 5. To exclude the possibility of any other 
solution we can use the following simple argument (valid for a = 2 l l a > 1, 
although a minor modification allows to deal with the case a < —1). 

Consider a function <f> G Coo(M) with cp(a) = /20(1) and || <p ||= 1, but 
arbitrary otherwise. Now for any positive integer we construct the function 
f {N) S Coo(M) defined for x G [a~ N ,a N ] by 

f( N \a n x) = fi n 0(x), x€[l,a], n = - N, -N + 1, . . . , N - 1, A, 
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and extended to I 



I so that || || 
C„(x)/W(x) dx 



1. Thus, 

,n+l 



Ux)fW(x) dx 



w G (l,a), 



-iV, 



(4.7) 



Then 



C^x)f( N \x)dx = (2N-l)I u 



and, if I w ^ 0, is an unbounded functional. Therefore /„, = 0, implying 
that the intersection of the support of (a with (0, oo) is empty. A similar 
argument rules out the negative semiaxis and finally one deduces that must 
be supported in {0}, i.e. it is the J- function that corresponds to [i = 1. 



As for the rest of the spectrum consider, for fi = e ie and a positive integer 



K, the following distribution E 



One can compute 



1 



c 



id log |x|/log |a| 



4 log if 



1 and 



if |z| #{1/K,K), 
if |z| G (l/-fC,if). 



log a 



(4.8) 



log if ' 

which shows that /i = e* e is in the spectrum of U . Therefore, the complex 
spectrum of (DT a )s = 2U is as stated in the second part of the proposition. 

As mentioned before, the only other possibility to obtain a fixed point 
with Re^4 = is having a = 1 and ip = ±ir/2. The fixed point in this case is 
5(x + e) with e = Im A The proof above can be applied to this case by simply 
shifting to — e the origin of the real line. □ 

We would like to discuss now the information about the behavior of the 
renormalization group transformation around the fixed points that can be ex- 
tracted from the study of the spectrum. We start by considering the fixed 
points corresponding to strictly stable Levy distributions. In this case the 
stable directions that lie on the domain of attraction of the fixed point cor- 
respond to | A s | < 1, i.e. s > a, while those with s < a give |A S | > 1 and 
therefore unstable perturbations of the fixed point. 

It is interesting to analyze the behavior of the characteristic function of 
the fixed point perturbed in the £s direction for small values of k. First note 
that 

p (k) = 1 - A\k\ a 6{k) - A\k\ a d(-k) + o(\k\ a ) 
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and that 



p {k) + e( s {k) 



1 - A\k\ a 9{k) - A\k\ a O(-k) + o(\k\ a ) for s > a, 
1 + B + k s 9(k) + B_\k\ s 9(-k) + o(\k\ s ) for s < a, 



which implies that the perturbations that do not change the behavior of the 
characteristic function around k = are precisely the stable directions. In 
this way we obtain an infinitesimal version of the well-known result |21j that 
the the p.d.f.s whose characteristic function behaves around k = as 



belong to the domain of attraction of the fixed point with characteristic func- 
tion 



We consider now the other fixed point, the Dirac delta function. From the 
results of Section 2 (see also below) one would expect to have an attractive 
fixed point at 5, at least for certain values of a. Therefore, it may seem sur- 
prising that the spectrum of the differential around the fixed point 5 has only 
eigenvalues with modulus greater than 1. The reason is that the eigenvectors 
with |A| < 1 are associated to the strong convergence of the iterated transfor- 
mation of the fixed point perturbed in that direction, while in the case of 5, 
strong convergence is impossible and we have convergence only in the weak 
topology. This fact can be shown by looking more closely at the behavior of 
the transformation around the fixed point. Consider a perturbation of 5 by 
h G C£o(K) with <|| h ||< 1. Then 

II T a (5 + h)-5 || = || 2h + h*h ||> 

> 2 || h || - || h* h ||> 2 || h || - || h \\ 2 >\\ h ||, (4.9) 

and all directions around the fixed point 5 are unstable in the strong topology. 
The situation in the weak topology is different and, as shown below, there are 
p.d.f.s in L 1 which under successive application of the renormalization group 
transformation converge weakly to the Dirac delta function. 

In order to complete the local picture obtained above and to substanti- 
ate previous statements, we shall discuss the behavior of the iteration of the 
transformation for different initial densities and different values of a. Let 



with < v < 2 and B ^ 0. Therefore, for fixed k and a = 2 l / a > 1 one has 



p{k) = 1 - A\k\ a 6(k) - A\k\ a 9(-k) + o{\k\ a ) 



p {k) = exp(-A\k\ a 6{k) - A\k\ a 9(-k)). 



p{k) = 1 - B\k\"8(k) - B\k\ u 9{-k) + o{\k\ v ) 



(4.10) 



T™p(k) = (1 - a- vn B\k\ u 9(k) - a- vn B\k\ u 9(-k) + o{ar vn )) 



2 
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and the limit of large n with k ^ reads 
lim T™p{k) = < 



1 for v > a, 

exp(-£|fc|"0(fc) - Elfcl^t-fc)) forzy = a, (4.11) 
for v < a. 



We therefore have the following proposition: 
Proposition 3: 

Take p € C^ )0 (R) with characteristic function given by \lj-10 ) and a = 2 1 / c 



Then, T™p converges weakly when n — > oo to 5 if v > a, to S a ,B if v = a and 
it does not converge if v < a. 

The result is a simple consequence of Levy's continuity theorem together with 
the behavior of the characteristic functions in (|4.1ip . □ 

With a similar computation we obtain for a = —2 l / a < — 1 the following: 

- If v > a, T a p converges weakly to the 6 function. 

- For v = a we have a limit two-cycle formed by the strictly stable densities 
S at B and S a The two members of the cycle are related by reflection 
with respect to the origin x i— > —x. 

- Finally, for v < a the sequence does not converge. 



In order to illustrate some of the advantages of our approach to the Gene- 
ralized Central Limit Theorem we will give now a simple proof, based on the 
above results, of the positivity of S a< A for a and A as in (|4.2p . The proof goes 
as follows: 

Take the following p.d.f. 

with q > 0, ci + C2 > 0, a G (0, 1). n is a positive real number, chosen so that 
|| p ||= 1. Standard estimates show that 

p(k) = 1 - A\k\ a 9{k) - A\k\ a 9(-k) + o(\k\ a ), 

where A = \A\e 1 ^ with 

\A\ = (c\ + c 2 2 + 2cic 2 cos(Tra)) 1 / 2 F(1 - ^ 

a 
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and 

tan0= tan( — ) with \<p\ < —. 

ci + C2 2 2 

On the other hand, from our previous results we have 

lim T^p = S at A, 

71— >00 

which, due to the fact that our transformation preserves positivity, implies 
that S a> A is positive. Note that taking different non-negative values for c\ and 
C2 we cover the whole range for A as described in (|4.2p , 

A slight modification of the initial p.d.f. and the estimate allows to deal 
with the remaining cases, i.e. a G [1,2]. 



5 Conclusions 

In this paper we have studied a simple instance of the renormalization group 
transformation on the space of probability densities. We have shown that 
by changing the scaling that governs the transformation one obtains different 
fixed points or limit two-cycles that correspond to the strictly stable laws of 
Levy and the limiting case of the Dirac 8 distribution. 

We have also studied the stability of the fixed points using the linear ap- 
proximation of the transformation around the fixed point. In this way we have 
derived a local version of the classical results about the domain of attraction 
of the stable laws. We have also shown that the stable or unstable character of 
the fixed point can be different if we consider the strong or the weak topology 
in the space of probability densities. 
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